A new phase structure in three-dimensional dynamical triangulation 
model 

T. Hotta'^, T. Izubuchi^ and J. Nishimura'^ 

'^Institute of Physics, University of Tokyo, Meguro-ku, Tokyo 153, Japan 
'^Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
"^Department of Physics, Nagoya University, Chikusa-ku, Nagoya 464-01, Japan 

We investigate a new phase structure of the three-dimensional dynamical triangulation model with an additional 
local term by Monte Carlo simulation. We find that the first order phase transition observed for the naive Einstein- 
Hilbert action terminates at a finite coefficient of the new term. The phase transition turns into second order at 
the endpoint, beyond which it becomes a crossover. 



1. Introduction 

The success of the dynamical triangulation 
model in two dimensions summarized by its 
equivalence to the Liouville theory in the con- 
tinuum limit motivated intensive studies of the 
model in higher dimensions through Monte Carlo 
simulation. It turned out, however, that with 
an action obtained by naive discretization of the 
Einstein-Hilbert action, the phase transition ob- 
served is of first order both in three dimensions 
and in four dimensions [|l|. In lattice theories in 
general, it is often preferable {e.g. improved ac- 
tions to reduce the finite lattice spacing effects) 
or even necessary {e.g. the Wilson term to elimi- 
nate doublers) to modify the lattice actions from 
that obtained by naive discretization of the ac- 
tion in the continuum. We may therefore hope 
to obtain a second order phase transition in the 
present case by adding a local term to the naive 
Einstein-Hilbert action. 

The new term we consider is motivated by a 
possible contribution from the measure for the 
path integral over the metric and has been first 
studied in four dimensions [|| . Recently this mod- 
ified action has been studied in three dimensions 
by the Monte Carlo renormalization group and 
a possibility was suggested that the first order 
phase transition terminates at a finite coefhcient 



of the new term ^. We clarify this possibility 
by high statistics Monte Carlo simulation using 
multi-canonical technique. 

2. The model 

The naive discretization of the Einstein-Hilbert 
action in the dynamical triangulation model in D 
dimensions gives Sq = —kqNo + koNjj, where 
iVo and N]j are the number of vertices and 
I?— simplices respectively in the simplicial man- 
ifold. We modify the action as 

s = So + ^lM, (1) 

M - ^ log [oH/(i^ + 1)] , (2) 

V 

where fi is the new coupling constant, and o{v) 
is the number of Z3— simplices that contain the 
node V. The summation in (|^) is taken over all 
the nodes in the manifold. This term can be 
considered as the contribution from the path inte- 
gral measure H^, (7(2;)"''/^, since o{v) can be inter- 
preted as the "local volume" , which corresponds 
to g{x) in the continuum. 

In the following we consider the canonical en- 
semble with fixed No given by the partition func- 
tion : 

Z{KQ,fi) = ^ exp{KoNQ ~ fiM) , (3) 

{T{No )} 

^Note that our fj, is the minus of the fj, in Ref. [pj 
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= ^z(iVo;^) exp(KoAfo), (4) 

No 

where the summation in (^) is taken over all pos- 
sible triangulations of Z3-dimensional sphere that 
contain Njj D-simplices. The z{Nq; /i) is the par- 
tition function for fixed Nq, which can be given 

by 

z(iVo;Ai)= J2 cxp(-/iM), (5) 

{T{Nd.No)} 

where the summation is now taken over all pos- 
sible triangulations that contain Nd i'-simplices 
and A'o vertices. The details of the numerical 
simulation will be reported elsewhere [Q. 

3. End point of the first order phase tran- 
sition 
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Figure 1: The microcanonical inverse tempera- 
ture K{No; /i) as a function of iVo for several /i's 
with the system size — 10, 000. 

In order to investigate the phase diagram with 
the modified action, we measure the microcanon- 
ical entropy s(Nq; ^jl) defined by 

s{No\^l) = log z{No-^i), (6) 

from which we calculate the microcanonical in- 
verse temperature K{No; fi) by 

K{No;^i)^s{No + ^;^i)-s{No;^i). (7) 

Precise measurement of the microcanonical en- 
tropy near the first order phase transition point 
is possible by using the multi-canonical technique. 



This was previously done for the unmodified ac- 
tion with /i = in Ref. and the critical 
exponents consistent with the first order phase 
transition have been successfully extracted by the 
finite-size scaling analysis. 

In Fig. |l|, we show the microcanonical in- 
verse temperature K{No;ii) for several values of 
/i. Note that the node distribution is given for 
each kq and ^ by z{No; fj,) exp^KoNo) and the dou- 
ble peak structure in the node distribution, which 
is a signal of the first order phase transition, can 
be probed directly by looking at the microcanon- 
ical inverse temperature K[Nq\ijl). The double 
peaks appear if and only if (i) K{No; n) has both 
local maximum and local minimum and (ii) kq 
is set to a value between the local minimum and 
the local maximum. One can see that the local 
minimum and the local maximum observed for 

= 1^ merge at ^* — 0.65, beyond which the 
K{No] fJ.) becomes a monotonous function of A'^o- 
This suggests that the first order phase transition 
ceases to exist at fj,*, beyond which it becomes a 
crossover. 

Free energy with \i term 
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Figure 2: The free energy /{Nq; ^, ko) as a func- 
tion of A'o near Kq for = and /i = /i* with the 
system size A^3 = 20, 000. 

For each fi below /z* , we define the critical cou- 
pling Kq of the phase transition as the position 
of the peak of the node susceptibility x(ko,m) = 
(A^q) — (Nq)'^ . The end point of the phase tran- 
sition fine is given by /i* = 0.65, Kq — 5.00 for 
N3 = 10, 000, and fi* = 0.80, n* = 5.515 for 
N3 = 20, 000. In Fig. |, we plot the free en- 
ergy f{No; ^j,;ko) = -s{No;^j.) ~ kqNq, which is 
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nothing but the minus log of the node distribu- 
tion, as a function of A^o near Kq for /it = and 
fi = fi* . The two local minima observed for fi = 
merge into one at n = /j*. The free energy curve 
at (kq, fi*) has a wide flat bottom, which suggests 
the existence of a massless mode. 

4. Fractal structure at the end point 

In order to examine a possible continuum limit 
at the end point, we measure the boundary 
area distribution p{A, r), which is the number of 
boundaries with the area A at the geodesic dis- 
tance r. The corresponding quantity in two di- 
mensions called the loop length distribution is cal- 
culated analytically and is found to possess a con- 
tinuum limit 1^. The scaling behavior expected 
from this result has been correctly reproduced by 
numerical simulation . 

In Fig. H, we plot the boundary area distri- 
bution at the end point of the first order phase 
transition line as a function of a; = A/r'^. 
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Figure 3: Boundary area distribution p{A, r) at 
the end point of the first order phase transition 
line with the system size — 20,000. 

One can see a reasonable scaling behavior. The 
power of r in the scaling variable x = A/r'^ im- 
plies that the fractal dimension is 3. We would 
like to remark that the scaling behavior we ob- 
serve is much better than the one that has been 
claimed to exist with the unmodified action S . 



5. Discussions 

Our results suggest that the first order phase 
transition becomes second order at the end point, 
where one can obtain a continuous theory. A nat- 
ural question to be asked is what the continuum 
theory is. One may suspect that we cannot con- 
struct any continuum theory in three-dimensional 
quantum gravity, since we have no physical de- 
grees of freedom. This argument is too naive, 
however. An example of field theory which has no 
physical degrees of freedom that is still completely 
well-defined is two-dimensional Yang-Mills the- 
ory 

The fact that we have to fine-tune two param- 
eters to obtain the continuum limit implies that 
there are two relevant operators around the fixed 
point. A possible interpretation of the continuum 
theory is therefore the gravity. The observa- 
tion that the fractal dimension at the fixed point 
is approximately three is also suggestive of this in- 
terpretation. The fractal dimension we extracted, 
however, is still preliminary, and the value might 
increase as we increase the system size. Needless 
to say, a lot of work is yet to be done to under- 
stand the nature of the continuum theory. 

In 4D, the situation is different from the 3D 
case in that singular vertices which are shared by 
huge number of Z?-simplices exist in the strong 
coupling phase |^ . We can add a local term in the 
action that eliminates these singular vertices, if 
one wishes, and then the situation becomes more 
or less the same . We therefore expect that the 
first order phase transition can be changed into 
second order by modifying the action. 
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